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ABSTRACT 


Experimental  data  on  the  crystal  properties,  second  virial  coefficients, 

and  viscosity  coefficients  of  Ne  , A , Kr  , Xe  , CH^  , N^  , CO  , 0£  , 

and  CO^  were  analyzed  for  the  purpose  of  obtaining  values  of  the  parameters  in 

g r (g  ot  ( i - r/r  ) / Y"m  *j 

the  exp-six  inter  molecular  potential,  j ~ G m — \^—J  J , 

For  gases  whose  molecules  are  spherical,  it  was  possible  to  reproduce,  with  a 
single  set  of  potential  parameters,  not  only  the  crystal,  second  virial,  and 
viscosity  data,  but  also  data  on  other  transport  properties  with  fair  accuracy. 
For  gases  whose  molecules  deviate  appreciably  from  spherical  symmetry  it  was 
necessary  to  choose  at  least  two  different  sets  of  potential  parameters  in  order 
to  reproduce  different  types  of  properties.  Such  behavior  was  taken  to  indicate 
the  inadequacy  of  the  assumptions,  made  in  the  fundamental  gas  theories,  that 
intermolecular  forces  are  central  and  that  inter  molecular  collisions  are  elastic. 
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Recent  calculations  of  the  transport  property  collision  integrals^  and 

; “r  j 

the  second  v trial  coefficients'  ' for  gases  whose  molecules  obey  ar  exp-six 
inter  molecular  potential  have  enabled  the  inter  molecular  potentials  of  helium 
ard  Lydrogou  to  be  .determined  from  experimental  values  of  the  viscosity  and 
second  virial  coefficients^.  , Th*  exp-six  potential  may  be  written  in  the  form 


q>(r) 


, ; i i } o 

6 U 


J . 


(1) 


where  Cj)  (r)  is  the  potential  energy  of  two  molecules  at  a separation  distance  r , 

£ is  the  depth  of  the  potential  energy  minimum,  rm  is  the  position  of  the 
minimum,  and  OC  is  a parameter  which  is  a measure  of  the  steepness  of  the 
repulsion  energy.  Corner^  has  shown  how  such  calculations  of  intern,  -lecular 
potentials  may  be  improved  by  the  use  of  crystal  data  in  addition  to  gas  properties 
data,  and  has  used  this  method  to  determine  the  potentials  of  neon  and  argon  from 
crystal,  second  virial  coefficient,  and  Joule -Thomson  coefficient  data.  In  the 
present  paper  we  extend  ths  application  of  Corner's  procedure  to  include  experi- 
mental transport  property  results  for  the  evaluation  of  the  intermolecular 
potentials  of  a number  of  simple,  non-polar  molecules. 


^ E.  A.  Maaon,  J,  Che  in.  Phvs.  22,  000  (1954). 

(2)  w.  E.  Rice  and  J.  O.  Hirschfelder,  J.  Chem.  Phys.  22,  000  (1954). 
^ E.  A.  Mason  and  W.  E.  -Rice,  J.  Chem.  Phys.  (to  be  published). 

(4)  ...  Corner,  Trans.  Faraday Soc.  44,  914  (1948). 
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DETERMINATION  OF  POTENTIAL.  PARAMETERS 


The  crystal  data  used  in  Corner's  method  are  the  lattice  distance  of  the 
crystal,  R , and  the  heat  of  sublimation,  Eq  , both  at  0°K.  The  heat  of 

o 

sublimation  is  corrected  for  the  zero-point  energy  by  taking  it  to  be  the  sum  of 
the  experimental  heat  of  sublimation  at  0°K  (in  units  of  energy /mole)  and  the 
quantity  (9NQk®jj/8)  , where  NQ  is  the  Avogadro  number,  k is  the  the 
Boltzmann  constant,  and  is  the  experimental  value  of  the  characteristic 

Debye  temperature  as  determined  at  low  temperatures.  Knowledge  of  the  two 
quantities  R and  Eq  permits  the  calculation  of  the  two  potential  parameters 
rm  and  € as  a function  of  the  third  parameter  o t by  solving  the  transcendental 
equation 


/ , *0-1?*) 
7 2(e-e)e 


.2 

■3  h ( oi-io) 
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?tr  me  ol 
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where  and  Cg  are  constants  tabulated  by  Lennard -Jones  and  Ingham^), 

( 4\ 

0 and  0'  are  functions  of  (ot  R*)  tabulated  by  Corner'***,  R*  is  the  reduced 
lattice  distance  &/rm  » h is  Planck’s  constant,  and  m is  the  mass  of  one 


(5) 


J.  E.  Lennard-Jones  and  A.  E.  Ingham,  Proc.  Roy.  Soc.  (London)  A107, 
636  (1925). 
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molecule.  The  value  of  the  parameter  € appearing  in  the  third  term  of  Eq.  (2) 
is  given  as  a function  of  o l and  R*  by  the  equation 


E = I ~ 72  & e 


(3) 


Eqs.  (2)  and  (3)  are  subject  to  the  conditions  that  the  crystal  is  cubic  and  that 
the  zero -point  energy  is  small^). 

The  procedure  for  finding  the  potential  parameters  rm  , & , and  od  is 
then  as  follows: 

(1)  For  a chosen  arbitrary  value  of  of,  , the  value  of  R*  which 

satisfies  Eq.  (2)  is  found  by  trial  and  error.  The  parameter  rm  is  found 

from  the  relation  r = R/R*  . 

m 

(2)  The  value  of  €■  is  next  calculated  directly  from  Eq.  (3). 

(3)  Another  value  of  oc  is  chosen  and  the  procedure  repeated. 

A pair  of  values  of  r and  6 is  thus  obtained  for  each  value  of  oc  selected. 

m 

(4)  The  value  of  oC  is  finally  found  from  experimental  values  of  the 
viscosity  and  second  virial  coefficients  by  the  method  of  translation  of  axes 
described  in  reference  (3).  In  this  method  a plot  of  experimental  quantities, 
for  example  log  B(T)  vs.  log  T where  B(T)  is  the  second  virial  coefficient 
and  T is  the  absolute  temperature,  is  moved  until  it  coincides  with  one  of 

a family  of  theoretical  curves;  the  amount  of  translation  parallel  to  each  axis 
determines  r^  and  € , and  the  particular  curve  chosen  from  the  family 


(6)  This  condition  prevents  the  use  of  crystal  data  in  determining  the  potential 
parameters  of  helium  and  hydrogen. 
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determine*  oC  . In  the  present  case  rm  and  & are  already  known  for  a given 
0 C , so  the  experimental  curve  is  moved  the  corresponding  amount  and  examined 
to  see  if  it  coincides  with  the  theoretical  curve  corresponding  to  the  same  ot  as 
do  r^  and  € . If  not,  a new  value  of  oc  is  selected  and  the  process  repeated 

until  the  values  of  o L coincide. 

The  above  procedure  obviously  over-specifies  the  potential  parameters; 
either  the  viscosity  or  the  second  virial  coefficient  would  be  sufficient  to 
determine  oC  , whereas  we  have  used  both.  This  fact  in  itself  affords  a test  of 
the  usefulness  of  the  potential  form  used,  since  an  unsuitable  potential  form 
would  be  unable  to  accomplish  the  simultaneous  fit  of  more  than  the  minimum 
number  of  properties  needed  to  evaluate  the  potential  parameters. 

Table  I lists  the  potential  parameters  found  by  the  above  procedure  for 
some  simple  non-polar  molecules,  and  Table  II  gives  a brief  summary  of  the 
crystal  data  used  to  determine  the  parameters.  For  comparison,  the  potential 
parameters  obtained  by  Corner'  ' from  crystal,  second  virial,  and  Joule -Thomson 
coefficient  data  for  neon  and  argon  for  the  Buckingham-Corner  potential^ 

(which  is  very  similar  to  the  exp-six  potential  except  at  small  values  of  the 
separation  distance)  are: 

^ R.  A.  Buckingham  and  J.  Corner,  Proc.  Roy.  Soc.  (London)  A189, 

118  (1947). 
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Table  I 

Exp-six  potential  parameters  for  some  simple,  non-polar  molecules, 


r 

Parameters 

a/ 

r 

m 

«/k 

(A) 

°K 

14.5 

3.  147 

38.0 

14.0 

3.866 

123.  2 

12.3 

4.056 

158.  3 

13.0 

4.450 

231.  2 

14.0 

4.206 

152.8 

16.2 

4.040 

113.  5 

4.011 

101.  2 

>17 

<4.  10 

>132 

17.0 

3.937 

119.  1 

R.  S.  Edwards,  Proc.  Roy.  Soc.  (London)  A119,  578  (1928);  M.  Trautz 
and  H.  E.  Binkele,  Ann.  Physik  5,  56 1 (1930):  M.  Trautz  and  R.  Zink, 
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TABLE  II 

Summary  of  crystal  data  used  for  the  determination  of  potential  parameters. 


Substance 

R 

(A) 

Eo 

(cal.  /mole) 

®D 

(°K) 

Ne 

3 :.o 

591 

64 

A. 

9 O 1 

M A 

TA  Oft 

LU  L7 

50 

Kr 

3.94 

2719 

63 

Xe 

4.  33 

3901 

55 

ch4 

4.  16 

2508 

78 

N2 

4.02 

1804 

68 

CO 

4.08 

2083 

79.5 

I 
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Ne; 

o t x 13.6, 

rm  = 3.  !6  A , 

e/k  = 37.  1 °K  ; 

A: 

o'  X 13.9, 

r = 3.87  A , 
m 

£/k  * 123.  2 °K 

The  agreement  with,  the  exp-six  parameters  o£  Table  I is  satisfactory. 

In  some  canes  it  was  not  possible  to  fit  all  measured  properties  witlT 
only  one  set  of  potential  parameters;  this,  failure  ic  discussed  in  the  next 
section.  It  was  necessary  to  extrapolate  the  tables  in  references  (1)  and  (2); 
the  supplementary  tables  obtained  by  the  extrapolations  are  given  in  the 
appendix. 


COMPARISON  WITH  EXPERIMENT 
A.  Spherical  Molecules 


Included  in  the  category  of  spherical  molecules  are  neon,  argon,  krypton, 
xenon,  and  methane.  Strictly  speaking,  the  theories  which  relate  the  second 
virial  coefficients  and  transport  properties  to  the  intermolecular  potential  are 
not  really  applicable  to  the  methane  molecule,  which  has  internal  degrees  of 
freedom.  However,  the  internal  degrees  of  freedom  do  not  seem  to  have  much 
effect  except  in  the  case  of  the  thermal  conductivity,  for  which  the  theory  fails. 

The  test  of  the  suitability  of  Eq.  (1)  as  a potential  form  is  its  ability  to 
reproduce  gas  properties  data,  since  the  crystal  data  are  automatically 
reproduced  by  the  procedure  used  to  determine  the  potential  parameters. 

Fig.  1 shows  a comparison  of  measured  values  of  the  second  virial  coefficient, 
B(T),  and  values  calculated  from  the  parameters  of  Table  I according  to  the 
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200  400 

T,  °K 


Comparison  of  observed  and  calculated  values  of  the  second  virial 
coefficients  of  spherical  molecules.  Solid  line:  exp-six  potential; 
broken  line:  Leonard -Jones  (12~6)  potential.  Observed  values: 

O neon,  A argon,  # krypton,  \7  xenon,  0 methane. 


equations 


B(T)  = b B*(ct  , T*)  , 
m 

b = (2  ft  13)  N r 3 , 
m ' ’ o m 
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(4) 

(5) 


where  B*(  ol  , T*)  is  a dimensionless  quantity  tabulated  in  reference  (2)  as 
a function  of  OC  and  T*  ( = kT / €.  ) , and  the  other  factors  are  as  previously 
defined.  Also  shown  are  the  values  of  B(T)  calculated  for  the  widely  used 
Lennar d- Jones  (12-6)  potential^ 


<((r) 


0 


(6) 


where  Cf  (r)  , 6 , and  rm  are  defined  in  connection  with  Eq.  (1).  Different 

values  of  the  parameters  rm  and  €.  of  the  L.  J.  (12-6)  potential  are  often 
available  for  the  same  substance,  depending  on  which  properties  were  used  to 
evaluate  the  parameters.  We  have  chosen  the  values  based,  at  least  in  part, 
on  crystal  data;  these  are  listed  in  Table  III.  The  agreement  between  the 
observed  and  calculated  values  of  B(T)  is  excellent  in  the  cases  of  neon  and 
argon  for  both  the  exp-  six  and  the  L.  J.  (12-6)  potentials.  The  agreement  is 
only  fair  for  methane  and  xenon;  for  xenon  the  exp-six  potential  fits  B(T)  better 
at  low  temperatures,  whereas  the  L.  J.  (12-6)  potential  fits  better  at  high 
temperatures.  The  agreement  in  the  case  of  krypton  is  good  for  the  exp-six 


(8)  All  tabulations  necessary  for  calculating  properties  of  gases  obeying  a 

' Lennard- Jones  (12--6)  potential  may  be  found  in  Hirschfelder,  Curtiss,  and 
Bird,  The  Molecular  Theory  of  Gases  and  Liquids  (John  Wiley  and  Sons, 
New  York,  1954). 


« * . wmrnvKim m t 1 mtm* 
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TABLE  NX 

Lennar d- Jones  (12-6)  potential  parameters  for  spherical,  non -polar  molecules. 


a.  MBestH  values  calculated  from  crystal  and  second  virial  coefficient  data 
by  J.  Corner  (reference  (4)  ). 

b.  Calculated  from  crystal  data  by  the  method  given  by  J.  Corner  in  Trans. 
Faraday  Soc.  35,  711  (1939). 
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potential,  bat  rather  poor  for  the  L.  J.  (12-6)  potential.  On  the  whole,  however, 
neither  potential  is  markedly  superior  to  the  other  as  judged  on  the  basis  of 
second  virial  coefficients. 

Fig.  2 shows  a similar  comparison  between  the  experimental  values  of 
the  viscosity  and  those  calculated  from  the  exp-six  parameters  of  Table  I 
according  to  the  equation 

j, 

i au.93  (Mr)* 

,0V~ 

m 


(3) 


(7) 


where  ^ is  the  viscosity  in  gm-cm- ^ -sec- ^ , M is  the  molecular  weight, 
rm  is  in  angstroms,  and  f^(o6,  T*)  and  XI  ' 2^*(  oc  , T*)  are  dimension- 
less functions  which  may  be  obtained  from  the  tabulations  in  reference  (1). 

Also  shown  are  the  values  calculated  from  the  L.  J.  (12  -6)  parameters  of 
Table  III.  All  the  experimental  values  were  corrected  to  a uniform  relative 
basis  in  which  the  viscosity  of  air  at  23°C  was  taken  as  1833.0  x 10"^  gm-cm"  * 
-sec"^  . Several  discrepancies  between  experiment  and  theory  are  apparent. 

The  measured  values  for  argon  deviate  badly  from  the  calculated  curves  at  high 
temperatures.  In  this  region  forces  of  repulsion  dominate,  and  the  slope  of  the 


experimental  curve  implies  an  exiremely  steep  repulsion  energy  at  close 


distances  of  approach  of  two  molecules.  However,  experiments  on  the  scattering 
of  high  velocity  argon  beams  in  argon  gas^  give  information  on  the  intermolecular 


(9) 


I.  Andur  and  E.  A.  Mason,  J.  Chem.  Phys.  (to  be  published). 
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Fit,,  2.  Comparison  of  ooservea  and  calculated  value*  of  the  viecosltle*  of 
spherical  molecule* , Solid  liae:  exp-eix  poteatial;  broken  line: 

Leaaard-Joaee  (12-6)  poteatial.  Obeerved  value*:  ©neon,  A argon, 

♦ krypton,  V xenon,  0 me  thane. 
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potential  in  this  region  of  strong  repulsion,  and  no  indication  of  such  a steep 
repulsion  is  found.  It,  therefore,  seems  likely  that  the  experimental  viscosity 
results  for  argon  are  low  at  high  temperatures.  The  deviations  at  low 
temperatures  for  argon  and  methane  are  greater  than  the  estimated  experimental 
error;  we  can  think  of  no  explanation  for  this.  The  rather  scanty  measurements 
of  the  viscosities  of  krypton  and  xenon  are  also  seen  to  deviate  from  the 
theoretical  curves.  Since  these  gases  are  difficult  to  obtain  pure,  it  is  possible 
that  some  light  gas  impurities  were  present.  Further  experiments  would  be 
helpful  in  clarifying  this  point.  Within  the  experimental  error,  the  exp-  six 
and  the  L.  J.  (12-6)  potentials  are  about  equally  successful  in  reproducing  the 
experimental  viscosities. 

A more  important  test  of  an  inter  molecular  potential  function  is  its 
ability  to  reproduce  properties  other  than  the  ones  used  in  determining  the 
potential  parameters.  Fig.  3 shows  a comparison  between  the  experimental 
values^®)  of  the  thermal  conductivity  and  those  calculated  from  the  exp-six 


(10) 


A.  Eucknn,  Fhysik.  Zeii.  12,  1101  (1911)  (A);  S.  Weber,  Ann.  Physik 
54,  437:  481  (l?17)(Ne  , A)  , Leiden  Comm.  Suppl.  42b  (1919)  (Ne); 

M.  Curie  and  A.  Lcpape,  Compt.  Rend.  19 3,  842  (1931),  J.  Phys.  et 
Rad.  2,  393  (1931)  (Ne,  A,  Kr,  Xe);  B.  G.  Dickens,  Pi ut . Roy.  Sou. 
(London)  Ai4jh  517  (1934)  (A);  W.  G.  Kannuluik  and  L.  H.  Martin, 
Proc.  Roy.  Soc.  (London)  A144,  496  (1934)  (Ne,  A);  W.  G.  Kannuluik 
and  E.  H.  Carman,  Proc.  Phys.  Soc.  (London)  B65,  701  (1952) 

(Ne,  A,  Kr,  Xe);  F.  G.  Knyes,  Project  Squid  Tech.  Mem.  MIT-  1 
(Oct.  1,  1952)  (A). 


Coiapari*oa  of  observed  aad  calculated  value*  of  the  thermal  coaductivitiee  o* 
•pherical  nolecule*.  Solid  line:  ;xp-*ix  potential;  broken  liae:  Leaaard-J' 
(12-6)  potential.  Observed  value*!  G aeon,  A argoa,  # krypton,  V xtm 
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parameters  of  Table  I according  to  the  equation 

1 

7 /m./  T s ) 

,0  * ' a * rv?) 


rt*  r 


/w 


(8) 


where  X is  the  thermal  conductivity  in  cal-cm"*-sec"*-deg“*  , f^i  o c , T*) 
ia  a dimensionless  function  tabulated  in  reference  (1),  and  the  other  terms  are 
as  defined  in  connection  with  Eq.  (7).  Also  shown  in  Fig.  3 are  the  values  of  X 
calculated  for  the  L.  J.  (12-6)  potential.  The  theory  from  which  Eq.  (8)  is 
derived  takes  no  account  of  energy  transported  in  the  molecular  internal  degrees 
of  freedom,  so  methane  is  excluded  from  this  comparison.  No  exact  method  for 
calculating  the  transport  of  energy  by  the  internal  degrees  of  freedom  has  yet 
boon  worked  out. 

Experimental  values  of  X are  much  more  difficult  to  determine  accurately 

than  values  of  7|  , but  within  the  experimental  error  the  exp-  six  and  L.  J.  (12-6) 

potentials  reproduce  the  measured  values  of  A about  equally  well.  The  agreement 

is  excellent  for  neon  and  argon,  only  fair  for  xenon,  and  poor  for  krypton.  It 

should  be  pointed  out  that  according  to  Eqs.  (7)  and  (8)  values  of  X and  for 

the  same  gas  are  not  independent,  but  instead  are  directly  proportional  (except 

/ < _ . 

for  the  functions  f and  f l5)  , which  are  very  nearly  equal).  Thus 
comparison  of  observed  and  calculated  values  of  X is  not  an  independent  test 
of  the  suitability  of  potential  parameters  derived  from  observed  values  of  , 
but  is  more  a check  on  the  consistency  of  the  measured  values  of  X and  ^ 
according  to  the  fundamental  theory  from  which  Eqs.  (7)  and  (8)  are  derived. 
Comparison  of  Figs.  2 and  3 shows  that  such  measurements  are  not  always 
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Fig.  4.  Comparison  observed  and  calculated  values  of  the  self-diffusion 

coefficients  of  spherical  molecules.  Solid  line:  exp-six  potential; 
broken  line:  Lennard-Jones  (12-6)  potential.  Observed  values: 

O neon,  A argon.  • krypton,  V xenon,  □ methane.  Curves 
a,  b,  c,  d,  and  e represent  neon,  argon,  krypton,  xenon,  and 
methane,  respectively. 
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calculated  for  the  exp -six  potential  and  the  L.  J.  (12-6)  potential.  The  agreement 

is  excellent  for  argon,  fair  for  neon,  and  rather  poor  for  krypton,  xenon,  and 

methane,  for  which  the  scatter  of  the  experimental  points  is  large.  Within  the 

scatter  of  the  observed  values,  both  potentials  are  zbeui  equally  satisfactory. 

A final  quantity  suitable  as  a test  of  potential  parameters  is  the  thermal 

diffusion  ratio,  . For  a binary  mixture  of  heavy  isotopes,  the  m-th 

approximation  to  the  thermal  diffusion  ratio,  k , may  be  written  as 

L T J m 


x,x 


12 


(10) 


correct  to  the  first  order  of  the  small  quantity  (M^  -M^/fMj+M^)  , where 
Mj  and  are  the  molecular  weights  of  the  two  isotopes,  and  Xj  and 

are  their  mole  fractions  in  the  mixture.  The  quantity  f k*  j m is  a reduced 
thermal  diffusion  ratio,  and  is  tabulated  in  reference  (1)  for  the  exp-six 

a 0e 

potential  a.3  a function  of  oi  and  T for  m = 2 . A similar  function,  designated 

as  R-p  , is  tabulated  for  the  L.  J.  (12-  6)  potential  to  a first  approximation  in 

j.  r * i 

reference  (8k  the  relation  between  krp  and  R js:  j k I - (10  5/1  18)R  . 

Experimental  values  of  have  been  obtained^  ^ for  neon,  argon,  and  methane 

by  measurement  of  the  thermal  separation  of  isotopic  mixtures  between  two 
absolute  temperatures,  T and  T'  . An  experimental  value  ol  k^,  is  ihus  a 


<12)  A.  O.  Nier,  Phys.  Rev.  56,  1009  (1939)  (CH4)  , Phys.  Rev.  57,  338  (1940) 
(Ne);  L.  G.  Stier,  Phys.  Rev.  62,  548  (1942)  (Ne,  A),  A.  K.  Mann,  Phys. 
Rev.  73,  413  (1948)  (Ne,  A);  A.  N.  Davenport  and  E.  R.  S.  Winter,  Trans. 
Faraday  Soc.  47,  1160  (1951)  (CH4). 
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itiiiitt  value  over  the  temperature  range  of  T to  T’  , but  It  may  be  shown^1^ 
that  such  a mean  value  is  equal  to  the  actual  value  at  an  intermediate  temperature, 
T , given  by 

tt» 

T = in  (T*/T)  . (11) 

T*  - T 

Fig.  5 is  a comparison  o£  the  experimental  and  calculated  values  of  the  reduced 

♦ — 

thermal  diffusion  ratio,  k^,  , as  a function  of  T . The  agreement  for  argon  and 
methane  is  slightly  better  for  the  exp-six  potential  than  for  the  L.  J.  (12-6) 
potential;  for  neon  the  agreement  is  poor  for  both  potentials.  The  continued 
increase  of  for  neon  at  high  temperatures  implies  a very  steep  repulsive 
potential^^;  the  existence  of  such  a steep  repulsion  is  supported  neither  by  the 
viscosity  data  shown  in  Fig.  2 nor  by  results  obtained  from  scattering  of  beams 
of  high  velocity  neon  atoms  by  neon  gas^.  Further  measurements  of  for 

neon  at  high  temperatures  would  seem  to  be  desirable. 


H.  Brown,  Phys.  Rev.  58,  661  (1940).  See  also  A.  N.  Davenport  and 
E.  R.  S.  Winter,  Trans.  Faraday  Soc.  47,  1160  (1951). 

T.  Kihara,  University  of  Wisconsin  Naval  Research  Laboratory  Report 
WIS  OCR - 7,  June  5,  1953  (to  be  published  in  Rev.  Mod.  Phys.). 
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Fig.  5*  Comparison  of  observed  and  calculated  values  of  tha  reduced  thermal  diffueiom  ratios 

for  spherical  molecules.  Solid  liae:  exp-six  potential;  broken  line:  Leaaard-Joaes 

(12-6)  potential.  Observed  value e:  O aeon,  A ergon,  □ methaae.  Curve*  a,  b, 

and  c represent  neon,  argon,  aad  methaae,  respectively. 
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Included  in  this  category  are  nitrogen,  carbon  monoxide,  oxygen,  and 
carbon  dioxide.  There  are  of  course  many  other  non-spherical,  non-polar 
molecules,  but  only  for  these  four  does  there  exist  a sufficient  quantity  of 
experimental  data  to  justify  the  assignment  of  three  disposable  parameters 
( & , rm  , G ) on  the  basis  of  the  measured  properties. 

The  L.  J.  (12-6)  potential  has  not  been  completely  successful  in  reproduc- 
ing the  measured  properties  of  gases  composed  of  non-spherical  molecules,  the 
principal  difficulty  being  that  different  potential  parameters  are  required  to  fit 
different  properties  of  the  same  substance.  The  question  arises  -whether  the 
fault  is  primarily  that  the  potential  form  is  inadequate,  or  that  the  basic  theory 
is  inapplicable  because  of  its  assumptions  of  spherically  symmetric  intermolecular 
potentials  and  elastic  intermolecular  collisions.  Some  light  might  be  thrown  on 
this  question  by  the  use  of  the  exp-six  potential,  which  is  more  realistic  and 
contains  more  disposable  parameters  than  the  L.  J.  (12-6)  potential.  If  the 
exp-six  potential  fails  for  non-spherical  molecules  in  the  same  way  as  does  the 
I*.  3m  (12-6),  it  might  reasonably  be  supposed  that  the  trouble  lies  primarily 
with  the  basic  theory  rather  than  with  the  potential  form. 

It  is  not  necessary  to  look  far  to  see  that  the  exp  six  potential  is  not 
entirely  successful  in  dealing  with  non-spherical  molecules.  Oxygen  and 
carbon  dioxide  cannot  be  handled  by  the  method  used  for  spherical  molecules 
because  they  do  not  crystallize  in  the  cubic  form.  Clearly  the  non-spherical 
nature  of  the  molecule  has  a strong  effect  on  the  crystal  properties,  but  it  is 
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not  even  possible  to  fit  both  the  second  virial  and  the  viscosity  coefficients  of 
oxygen  and  carbon  dioxide  with  a single  set  of  potential  parameters.  The  best 
that  can  be  done  is  to  assign  different  parameters  for  each  property,  thereby 
reducing  the  exp-six  potential  to  the  status  of  an  unnecessarily  complicated 
interpolation  function. 

The  situation  is  somewhat  better  for  nitrogen  and  carbon  monoxide,  which 
crystallize  in  the  cubic  system.  Now  it  is  possible  to  fit  both  the  crystal 
properties  and  the  second  virial  coefficients  with  a single  set  of  potential 
parameters,  but  impossible  to  fit  both  the  crystal  properties  and  the  viscosities 
for  any  reasonable  value  of  the  parameter  oi.  . This  behavior  is  shown  by  our 
assignment  of  two  sets  of  potential  parameters  for  nitrogen  and  carbon  monoxide 
in  Table  I.  The  situation  is  more  clearly  illustrated  in  Figs.  6 and  7,  where 
the  experimental  second  virial  and  viscosity  coelficients  of  nitrogen  are  compared 
with  the  values  calculated  from  the  two  sets  of  parameters  m Table  I.  Each  set 
of  parameters  accurately  reproduces  the  experimental  results  for  which  the 
parameters  were  evaluated,  but  fails  to  fit  the  other  experimental  results. 

Although  the  exp-  six  parameters  derived  from  viscosity  do  not  reproduce 
the  crystal  properties  and  second  virial  coefficients  for  non-3pher ical  molecules, 
they  may  be  fairly  successful  in  reproducing  other  transport  properties,  as  is 
illustrated  in  Fig*.  8 and  9,  where  experimental  self -diffusion  coefficients^ ^ 
and  thermal  diffusion  ratios  of  nitrogen  are  compared  with  the  values 

(15)  E.  B.  Win n,  Phys.  Rev.  74,  698  (1948);  80,  1024  (1950);  E.  R.  S.  Winter, 

Trans.  Faraday  Soc.  47,  342  (1951). 

A.  K.  Mann,  Phys.  Rev.  73,  413  (1948);  A.  N.  Davenport  and  E.  R.  S. 

Winter,  Trans.  Faraday  Soc,  47,  1160  (1951). 
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Comparison  of  thu  observed  values  of  the  second  virial  coefficient  of 
nitrogen  with  the  -value*  calculated  for  the  exp-six  potential.  Solid  line; 
calculated  using  the  parameters  derived  from  crystal  properties  and 
second  virial  coefficients;  broken  line:  calculated  using  the  parameters 
derived  from  viscosity. 
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Fig.  7.  Comparison  of  the  observed  viscosities  of  nitrogen  with  those  calculated 
for  the  exp-six  potential.  Solid  line:  calculated  using  the  parameters 

derived  from  crystal  properties  and  second  virial  coefficients;  broken 
lme:  calculated  using  the  parameters  derived  from  viscosity. 


Fig.  8.  Comparison  of  the  observed  self-diffusion  coefficients  of  nitrogen  with 
those  calculated  for  the  exp -six  potential.  Solid  line:  calculated  using 
the  parameters  derived  from  crystal  properties  and  second  virial 
coefficients;  broken  line:  calculated  using  the  parameters  derived  from 
viscosity. 
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Fig.  9.  Comparison  of  the  observed  reduced  thermal  diffusion  ratios  o 
nitrogen  v/ith  those  calculated  for  the  exp-six  potential.  Solid 
line:  calculated  using  the  parameters  derived  from  crystal 
properties  and  second  virial  coefficients ; broken  line:  calculated 
using  the  parameters  derived  from  viscosity. 
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calculated  from  the  two  sets  of  parameters  for  nitrogen.  The  parameters 
derived  from  the  viscosity  give  better  agreement  than  those  derived  from 
crystal  and  second  virial  coefficient  data.  Thus  if  it  is  necessary  to  calculate 
some  property  of  a gas  composed  of  non-spherical  molecules,  it  is  probably 
best  to  use  parameters  derived  from  equilibrium  measurements  if  the 
desired  property  is  am  equilibrium  one,  and  to  use  parameters  derived  from  the 
viscosity  for  calculating  transport  properties. 

SUMMARY 

The  three-parameter  exp  six  inter  molecular  potential  has  bsen  show-- 

to  be  able  to  predict,  with  a single  set  of  potential  parameters,  the  properties  of 

gases  composed  of  spherical  molecules  with  fair  accuracy.  For  the  gases 

considered  here,  however,  the  simpler  two  parameter  L,.  J.  (12-6)  potential 

f 3) 

is  scarcely  inferior  (in  contrast  to  the  case  of  helium  and  hydrogen'  ').  For 
gases  whose  molecules  deviate  appreciably  from  spherical  symmetry,  a single 
set  of  potential  parameters  is  not  able  to  reproduce  all  the  measured  properties. 
Since  the  form  of  the  exp--six  potential  is  physically  realistic  and  quite  flexible, 
such  behavior  would  seem  to  indicate  that  the  assumptions  of  central  inter - 
molecular  forces  and  elastic  intermolecular  collisions  made  in  the  basic 
theory  are  inadequate  to  describe  the  behavior  of  most  real  gases. 
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APPENDIX 

The  tabulations  of  transport  property  collision  integrals  and  second 
virial  coefficients  given  in  references  (1)  and  (2)  for  the  exp-six  potential 
extend  from  od  = 12  to  15.  This  range  of  the  parameter  oc  is  sufficient  for  all 
spherical  molecules  so  far  investigated,  but  non-spherical  molecules  seem  to 
require  larger  values  of  OC  , We  have,  therefore,  extended  the  tabulations  by 
extrapolation,  with  a consequent  reduction  in  accuracy.  Whereas  the  tabulations 
for  CL  = 12  to  15  were  for  the  most  part  accurate  to  within  a few  tenths  of  one 
percent,  the  extrapolations  may  well  involve  uncertainties  amounting  to  several 
percent. 

The  reduced  second  virial  coefficients,  B*(  oc  , T*)  , defined  by  Eq.  (4) 
were  extrapolated  graphically  to  od  = 17.  The  values  of  B*(  oc  , T*)  are  nearly 
linear  in  (log  ) , so  that  the  extrapolation  is  simple.  Values  of  B*(  oc  , T*) 
are  given  in  Table  IV  for  06  = 16,  17. 

Some  of  the  reduced  collision  integrals,  X2  ^ ( oc  , T*)  , two  of 

which  appear  explicitly  in  Eqs.  (7)-(9),  were  extrapolated  to  oc  =.  16,  17. 

These  functions  are  nearly  linear  in  cL  . The  functions  X),  ^ ^ are  not 

tabulated  directly,  but  instead  the  related  functions, 

ZU'n|=  [t.(1  .6/0,  )]*Q  a ■ “>*.  <U> 


which  are  more  suitable  for  interpolation  purposes,  are  given  in  Table  V. 
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TABLE  IV 

The  reduced  second  virial  coefficients. 


B*(£*  , T*) 


T* 

CL  = 16 

= 17 

.40 

-8.60 

-8.  28 

.45 

-6.  68 

-6.44 

.50 

-5.40 

-5.  19 

.60 

-3.82 

-3.67 

.70 

-2.88 

-2.  77 

.80 

-2.  26 

-2.  16 

.90 

-1.  83 

-1.75 

1.0 

-1.  505 

-1.429 

1.2 

-1.063 

-1.000 

1.4 

-0. 772 

-0.717 

1.6 

-0.568 

-0.518 

1.8 

-0.420 

-0.375 

2.0 

-0. 302 

-0.259 

2.  5 

-0.  107 

-0.070 

3.0 

0.016 

0.048 

3.5 

0.099 

0. 128 

4.0 

0.159 

0.  188 

4.5 

0.  201 

0.  227 

5. 

0.236 

0.  261 

6 

0.  284 

0.  307 

7 

0.316 

0.  338 

3 

0.  338 

0.360 

9 

0.  354 

0.  375 

10 

0.366 

0.387 

20 

0.  397 

0.417 

30 

0.393 

0.412 

40 

0. 383 

0,402 

50 

0.373 

0.  393 

60 

0.  364 

0.  384 

70 

0.  356 

0.  376 

80 

0.348 

0.  368 

90 

0.  341 

0.  361 

100 

0.334 

0.  355 
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(lt  2)  - (2,  2 ) „(2,  3) 


- 17 


,(1,  1)  (l,  2)  „ (2,  2) 


1.  195 
1.  197 
1.  204 
1.  207 
l.  198 
1.  178 
1.  154 
1.  130 
1.  108 
1.088 
1.072 
1.048 
1.032 
1.023 
1.019 
1.017 
1.023 
1.036 
1.052 
1.069 
1.  103 
1.  135 
1.  165 
1.  193 
1. 219 
1. 242 
1. 284 
1.321 

t 1 C.A 

*•  „ V 71 

1.  384 
1.411 
1 . 470 
1. 520 
1 . 564 
1.60  3 
1 . 639 
1.671 
1. 730 
1.782 
1.828 
1.871 
1.910 


1.095 
1.099 
1.  107 
1.10  5 

1.082 
1.051 
1.021 
0.995 
0.974 
0.958 
0.946 
0.931 
0.925 
0.925 
0.928 
0.933 
0.951 
0.972 

0. 993 
1.014 
1.053 
1.088 

1.  119 
1.  147 
1.  172 
1.  195 
l.  236 
1. 272 

I 3A  A 

X • JU  “I 

1.  332 
1. 358 
1.415 
1.463 
1.506 
1. 544 
1. 578 
1. 610 
1.  668 
1. 719 
1. 765 
1.807 
1 . 847 


. 187 
. 160 
. 148 
. 125 
.099 
.074 
.050 
.029 
.012 
.998 
.987 
0.971 
0.962 
0.957 
0.956 
0.957 
0.966 
0.979 

0. 994 

1.010 

1.041 
1.070 
1.097 

1.  122 
1.  145 
1.  167 
1 . 20  5 
1.  239 


1,  £>IU 

l.  298 
1.  323 
1.  379 
1.427 
1.468 
1.506 
1. 540 
1. 571 
1. 627 
1.  677 
1.722 
1.763 
1.801 


1.055 
1.025 
1.010 
0.977 
0.943 
0.912 
0.888 
0.871 
0.858 
0.849 
0.844 
0.  839 
0.  840 
0.844 
0.850 
0.858 
0,  878 
0.899 
0.920 
0.940 

0. 976 

1. C07 
1.036 
1.061 
1.085 
1.  106 
1.  144 
1.  177 
1 . 207 
1.  233 
1. 258 
1.311 
1.357 
1. 397 
1 434 
1. 466 
1.497 
1 . 552 
1.601 
1.645 
1. 685 
1.723 


. 195 
. 190 
. 196 
.201 
. 196 
. 180 
1.  159 
1.  138 
1.  119 
1.  101 
1.087 
1.065 
1.052 
1.044 

1.041 

1.041 
1.050 
1.065 

1.083 
1.  102 
1.  139 
1.  175 
1.  208 
1. 238 
1. 266 
1.  292 
1.  3 39 
1.  379 
1.415 
1 448 
1.478 
1.  543 
1 . 598 
1.  b47 
1.689 
l.  728 
1 » 764 
1.828 
1.885 
1.936 
1.983 
2.026 


r(2»  3) 


1.095 
1.092 
1.  100 
i . iu2 

1.084 

1.058 

1.031 

1.007 

0.988 

0.974 

0.963 

0.950 

0.946 

0.947 

0.951 

0.958 
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In  Table  VI  are  given  the  values  of  the  following  ratios  of  collision 
integrals  which  occur  frequently: 

a,  = a (2.  2)*  / SI  (1’  *>*  , (13) 

c*  = xi  <*•  2>*  / n u.  d*  . 

Table  VI  also  gives  the  values  of  the  first  approximation  to  the  reduced  thermal 
diffusion  ratio,  ["  kr£  ] j , as  calculated  from  a simple,  accurate  formula 
derived' by  Kihara^^: 


[4]  . = « issii J!  . 

L 1J  1 16  A* 


(14) 


"This  expression  gives  values  of  k*  which  differ  from  the  values  calculated  from 
the  elaborate  expressions  for  | k*  j ^ in  reference  (1)  by  no  more  than  the 
uncertainties  involved  in  the  extrapolations  of  the  XI  ^ ^ . 

The  quantities  f(m)(a,  T*)  which  appear  in  Eqs.  (7) -(9)  are  functions 
of  the  reduced  collision  integrals,  £1  ( ^ , and  are  used  for  calculating 

higher  approximations  to  the  transport  coefficients.  The  superscript  indicates 
the  order  of  approximation  used  in  computing  these  functions.  Kihara^^  has 
given  simple  expressions  for  these  functions  which  differ  from  the  more 
elaborate  expressions  used  in  reference  (1)  by  less  than  the  uncertainty  involved 
in  the  extrapolations  of  the  Xi  (•*■  • . The  values  of  the  f s were  therefore 

calculated  according  to  Kihara's  expressions: 


AZ)  . 1 f 

tv,  = 1 + 


-!]\ 


\ ~ a4**”* 

(S,3iA 

f(2)  ...  J Z r __  . 7]" 

* a*-*9*  *]  ’ 


(15) 


e(2)  = , + 


and  are  tabulated  in  Table  VII. 

(17)  •j>>  Kihara,  Imp  erf  e ct  Gases,  originally  published  in  Japanese  (Asakusa 
Bookstore,  Tokyo,  1949)  and  translated  into  English  by  the  I).  S.  Office 
of  Air  Research.  Wright-Patterson  Air  Force  Base. 
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The  functions  A* 


TABLE  VI 

r * 7 

and  the  reduced  thermal  diffusion  ratio.  I k Z I . . 


Ot  = 17 


1.093 


1.  118 


0. 889 
0.886 
0.880 
0.  866 
0.852 
0.844 
0.840 
0.840 
0.842 
0.845 
0.849 
0.859 
0.  868 
0.877 
0.885 
0.891 
0.905 
0.914 
0.921 
0.927 
0.934 
0.938 
0.940 
0.942 
0.944 
0.944 
0.945 
0 . 946 
0.946 
0.946 


lkT  3 1 


0.311 
0.286 
0. 251 
0.  170 
0.097 
0.052 
0.0  34 
0.0  32 
0.042 
0.060 
0.081 
0.  130 
0.  179 
0.223 
0.263 
0.299 
0.368 
0.417 
0.453 
0.478 
0.513 
0.531 
0.  544 
0.  550 

n e c 7 

\j  • ~j  j i 

0.560 
0.  564 
0.  56  3 
0.564 
0. 564 
0.564 
0.563 
u.  563 
0.564 
0.  564 
0. 564 
0. 567 
0.569 
0.  572 
0.574 
0.578 
0. 582 


.006 
.036 
.050 
.070 
.088 
.099 
. 104 
. 106 
. 105 
. 104 
. 102 
.097 
.093 
.091 
.089 
.088 
.087 
.087 
.089 
1.091 
1.095 
1. 098 
1.  101 
1.  104 

1.  j.06 
1.  108 
1.  Ill 
1.  113 
1.  114 
1.  116 
1.  117 
1.  119 
1.  120 
1.  121 
1.  121 
1.  122 
1.  123 
1.  124 
1.  124 
1.  124 
1.  125 
1.  125 


0.953 

0.954 

0.954 

0.955 

0.956 

0.957 


k) , 


0.311 
0.287 
0.262 
0. 193 
0.  127 
0.084 
0.066 
0.064 
0.074 
0.091 
0.111 
0.  158 
0.206 
0.249 
0.289 
0.323 
0.  391 
0.439 
0.475 
0.500 
0.  534 
0.  553 
0.566 
0.572 
0.580 
0.  582 
0.  588 
0.587 
0.590 
0.589 
0.591 
0.591 
0.  591 
0.593 
0.595 
0.596 
0.590 
0.602 
0.605 
0. 609 
0.613 
0.618 
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TABLE  VIJ 

Functions  for  calculating  higher  approximation^  to  the,  transport  coefficients. 


oL  ~ 16 

Oi  = 17 

T* 

f'2) 

\ 

f(Z) 

D 

£«’  23 

f(2) 

1 

f(2) 

D 

o 

1 . 00 17 

1.0026 

1.0020 

1. 0017 

1,0026 

1.0020 

. 1 

1 . 00  1 8 

1.0028 

1,0018 

1. 0018 

1.0028 

1.0018 

■} 

1.0020 

1.0031 

1.0014 

1. 0020 

1.0030 

1.0015 

.3 

1.0016 

1.0025 

1.0007 

1. 0017 

1.0027 

1. 0009 

.4 

1.0008 

1.00  12 

1.000  2 

1.0010 

1.0016 

1.0004 

.5 

1.0003 

1.0005 

1.0001 

1. 0005 

1. 0007 

1.0002 

.6 

1.  0001 

1.0001 

1.0000 

1.0002 

1.  0003 

1.0001 

.7 

1.0000 

1.0000 

1.0000 

1.0001 

1.0002 

1.0001 

.8 

1.0000 

1.0000 

1.0000 

1. 0001 

1.0001 

1.0001 

.9 

1.0000 

1. 0000 

1. 0001 

1.0001 

1.0001 

1.0002 

1.0 

1.0001 

1.0001 

1.0002 

1.0001 

1.0002 

1.0003 

1.2 

1.0002 

1.000  3 

1.0004 

1.0003 

1.0005 

1.0006 

1.4 

1.0004 

1.0007 

1.0008 

l. 0006 

1.0009 

1.0010 

1.6 

1.0008 

1.0013 

1.0012 

1. 0010 

1.0016 

1.0015 

1.8 

1.0013 

1.0019 

1.0016 

1. 00 

1.0023 

1.0020 

2.6 

1.0017 

1.0027 

1.0021 

1.0020 

1.0030 

1.0024 

2.5 

1.0029 

1.0045 

1. 0032 

1.0032 

1 .'0049 

1.0036 

3.0 

1.0039 

1.0060 

1. 0041 

1.0042 

1.0065 

1.0045 

3.  5 

1 0047 

1.0073 

1.0049 

1.0051 

1.0079 

1.00  53 

4 

1.0054 

1.0083 

1 .00  54 

1.0057 

1.0089 

1.0059 

5 

1.0062 

1. 0096 

i .0063 

1. 0066 

1.0103 

1.0068 

6 

1.0067 

1.0  105 

1.0068 

1. 0072 

1.0111 

1.0013 

7 

1.0071 

1.01  10 

1.00  72 

1.0075 

1.0117 

1. 0077 

8 

1.  0071 

1.0113 

1.0073 

1.0077 

1.0120 

1.0079 

9 

1.00  74 

1.0115 

1 .0076 

1.0078 

1.0122 

1.0081 

10 

1. 0074 

1,0116 

1.0077 

1.00  79 

1.0123 

1.0082 

12 

1.0075 

1.0117 

1 . 0078 

1.00  80 

1.0124 

1.C084 

14 

1.0075 

1.0117 

1. 00  78 

1.0080 

1.0124 

1.0084 

16 

1. 0075 

1.0117 

1.00  79 

1. 00  80 

1.0124 

1.0085 

18 

1. 0075 

1.0117 

1.0079 

1.0080 

1.0124 

1.0085 

20 

1.0075 

1.0117 

1 0079 

1.0079 

1.0123 

1.0086 

25 

1. 0075 

1.0117 

1.0079 

1.0079 

1.0123 

1.0086 

30 

1. 0075 

1.0117 

1.0079 

1.00  79 

1.0123 

1.0086 

35 

1.0076 

1.0117 

1.0080 

1. 0079 

1.0123 

1.0087 

40 

1.0076 

1.01  J6 

1.0080 

1. 0080 

1.0124 

1.0087 

45 

1.0076 

1.0119 

1.0080 

1.00  80 

1.0124 

1.0088 

50 

1. 0077 

1.01  19 

1.0081 

1.0080 

1.0124 

1.0088 

60 

1.0078 

1.0121 

1.0082 

1.0081 

1.0126 

1.0090 

70 

1.0079 

1.0123 

1.0083 

1.0082 

1.0128 

1.0091 

80 

1.0081 

1.0125 

1.0084 

1.0083 

1.0130 

1.0Q92 

90 

1.0031 

1.0126 

1.0085 

1.0084 

1.G131 

1.0093 

100 

1.0082 

1.0128 

1.0087 

1.0085 

1.0133 

1.0095 

